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I Abstract 

The instability of non-homoentropic axisymmetric flow of perfect fluid with respect to non-axisymmetric in- 
, finitesimal perturbations was investigated by numerical integration of hydrodynamical differential equations in two- 
' dimensional approximation. The non-trivial influence of entropy gradient on unstable sound and surface gravity 
waves was revealed. In particular, both decrease and growth of entropy against the direction of effective gravita- 
tional acceleration g^f / give rise to growing surface gravity modes which are stable with the same parameters in the 
case of homoentropic flow. At the same time increment of sound modes either grows monotonically while the rate 
, of entropy decrease against g^ff gets higher or vanishes at some values of positive and negative entropy gradient 
CIh in the basic flow. The calculations have showed also that growing internal gravity modes appear only in the flow 
i unstable to axisymmetric perturbations. At last, the analysis of boundary problem with free boundaries uncovered 
, that's incorrect to set the entropy distribution according to polytropic law with polytropic index different from 
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adiabatic value, since in this case perturbations don't satisfy the free boundary conditions. 
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1 Introdution 



This paper is natural extension of the two previous, in which the dynamical instability of axisym- 
metric flows was considered in approximation of incompressible perfect fluid (Zhuravlev & Shakura, 
2007a - hereafter ZSl) and with finite compressibility (Zhuravlev & Shakura, 2007b - hereafter 
ZS2). Papaloisou and Pringle (1984, 1985, 1987) were the first who investigated the instability to 
non-axisymmetrical perturbations being under discussion in this paper. Their discovery entailed 
numerous explorations, which are mentioned partially in the review by Papaloisou & Lin (1995). 
The growth of perturbations in axisymmetric flows with free boundaries is an important problem 
for astrophysics and is of special interest in the view of question how angular momentum transports 
outwards in accretion flows. Though the realistic astrophysical conditions force us to include the 
effects of radiative heating and cooling, magnetic fields, relativity and so forth, it's still necessary to 
study the accretion flow as a shear flow in the simple hydrodynamical approximation. In the most of 
such investigations the angular velocity profile was set by a power law fl oc r^'^ with 1.5 < g < 2.0, 
which means that Qe/f attains the greatest value at the boundaries. In ZSl and ZS2 alone with the 
power law we considered another fl profile which implies that geff = at the boundaries. It was 
revealed that the absence of effective gravity at the edge of the flow essentially affects the growing 
perturbation modes. In particular, it causes the disappearance of unstable sound modes so that the 
flow with angular velocity proflle being close to the Keplerian one becomes stable to infinitesimal 
perturbations. Additionally, the special attention was paid to the infiuence of vorticity gradient 
on the instability. In the present paper we are going to study how stratification of the basic fiow 
affects the growing sound and surface gravity waves and make an attempt to find unstable internal 
gravity waves. Following the other authors we consider the barotropic basic fiow, i.e. the entropy 
distribution is set in such way that isopycnic surfaces coincide with isobaric ones. It allows to regard 
angular velocity as function of radial direction. Discussing the instability of stratified axisymmetric 
flows we should refer, for instance, to the paper by Prank & Robertson (1988) who considered the 
instability of tori with random initial small perturbations and Kojima et al. (1989) who studied 
toroidal as well as cylindrical flows. The latter are two-dimensional analogue of real flows in point- 
mass gravitational potential. We should note here that in both cases the obtained results proved to 
be similar and 3-dimensional growing modes of perturbations in the toroidal flow turned out to have 
no essential dependence on vertical direction. The authors explained this by the fact that Reinolds 
stress didn't contain vertical perturbation velocity for barotropic conflgurations with Q(r). Further, 
Glatzel (1990) considered the instability of cylindrical and plane-parallel flows in approximation of 
small extent of a shear layer. In order to exclude the growing sound and surface gravitational modes, 
he assumed an incompressible fluid and rigid boundaries so that stratiflcation was entered by the 
density proflle. The instability he found was explained as the result of over-reflection and coupling 
of internal gravitational modes. Later Ghosh & Abramowicz (1991) studied the cylindrical flow con- 
sisting of two fluids with different constant densities. The configuration was Rayleigh- Taylor stable. 
Alone with modified branch of growing surface gravity mode being consequence of the free bound- 
aries (Blaes & Glatzel, 1986) they found another branch of increments appeared owing to density 
discontinuity. This instability is caused by growing internal gravitational mode, which, however, is 
fully analogous to the growing surface mode, since the boundary between two fiuids differs from the 
free boundary just by the finite ratio of densities. 

We should also mention about the results of Lovelace et al. (1999) and Li et al. (2000) where in the 
2-dimensional geometry the instability of thin Keplerian disks with the local maximum of entropy was 
investigated. In particular, the local dispersion relation for non-axisymmetric perturbations similar 
to Rossby waves dispersion relation was derived. At last, Klahr & Bodenheimer (2003) studied the 
instability of Keplerian disks with entropy decreasing outwards. 

As it's well known, the stability of stratified flows is governed by Richardson criterium. Initially it 
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was formulated for plane-parallel flows (see, for instance, Howard, 1961). In this case, the sufficient 
condition for stability is Ri > 1/4 in each point, where Ri - is the Richardson number. The gen- 
eralization of the Richardson criterium for baroclinic rotating flows, when angular velocity depends 
also on the vertical direction, was found by Fujimoto (1987) for incompressible fluid and by Hanawa 
(1987) for compressible fluid. 

In the present work we follow other authors and consider the simple 2-dimensional geometry, i.e. 
assume the cylindrical flow and perturbations with zero vertical velocity. In the flrst section we'll 
write down the basic equation for inflnitesimal perturbations, formulate the free boundary conditions 
and set the main features of the basic flow - angular velocity proflle and entropy distribution. Then 
we'll touch on a numerical method, in particular, discuss the question when the perturbations can 
satisfy the boundary conditions and flnally will present our results. 

2 Equation and boundary conditions. 

2.1 Derivation of main equation 

An isentropic motion of fluid is governed by the system of equations: (see, for example, Brekhovskikh 
& Goncharov, 1982): 

, Vp „^ 

+ vV V = - V$, 

ot p 

% + V(pv) = 0, (1) 



dt [p\ 



0, 



where $ is an external gravitational potential; The last equation represents the conservation of 
entropy in the particle of fluid with the perfect gas equation of state - p = K e^^^^'p"', with 7 = Cp/cy 
- adiabatic index, S* c„ - speciflc entropy and heat capacity for constant volume. 
The linearization procedure gives the equations for small perturbations: 

95v , / r V(5p 5p Vp 

— - + vV (5v + 5vV V = + — — , 

at P P P 

^ + V(pH + V(5pv) = (2) 



^f6p_ 6pp\ 
dt\p^ P p^ ^ ' 



f^-7^^U(^vV)i^ = 0, 

p py p^ 



where the third one implies the relation between perturbations of density and pressure. 

In our approach all quantities are functions of r, so for Eulerian perturbations with 6vz = in 
the cylindrical coordinates the system ([2]) will take the form: 

dSvr ^ n ^^^^ 2CJv - ^ ( ^pdp\ 
dt dip p\dr p dr j ^ 

-^ + n^ + -—{nr^)6vr = 3 

dt d(p r dr rp dtp 

d6p ^d6p 1 d . ^ . 1 d . ^ . 
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We'll be looking for a solution in the form of normal modes: 

here m - is the azimuthal wavenumber, uj = uj^. + iui - complex frequency, which for growing modes 
has the positive imaginary part, i.e. increment of mode. Possible uj has to be determined from the 
solution of the boundary problem coming to the integration of ordinary differential equation with 
the boundary conditions in the boundary points ri and r2 for functions v, p p, which, in general, 
are complex. 

The substitution (jlj) into ([3]) yields: 

^ • , • n\- on- ^ ( P^P\ 

i—tUi} + imiljVr — iiifip = ;- , 

p \dr pdr J 

{—tuj + tm\L)v^ -\ ;- I Vr = p (5) 

r dr ^ ^ rp 

1 (i _ imp _ 

[—luj + im\l)p H rvP'^r) H = 0, 

r dr r 

(—iu + imVt) \p — 7P- I + p^Vr— I — I = 
V P) dr \p^J 

From the last expression in ([5]) we'll get the relationship between p and p: 

P ? 

- = —--ir , (6) 

P « P 9eff 

where a? = •yp/ p - is the squared sound velocity, 

2 GM I dp 
9eff = i^ --^)^ = -Tr 

- the effective gravitational acceleration and = ivr/ {u> — rriQ) gives the radial dependence of radial 
component of the Lagrangian displacement : ^^(r, t) = ^r(r) e^'^'^*-'"^). At last, 

^ [-pTr-^) 

- is the squared buoyancy, or Brunt- Vaasala, frequency which is the characteristic internal oscillation 
frequency in medium at rest due to the entropy gradient and gravitational force. Clearly, for constant 
entropy dp = a'^dp, N'^ = 0, and the relation between the Eulerian perturbations of density and 
pressure comes to the expression one has for homoenropic configuration, when the perturbed flow as 
well as the unperturbed one are barotropic: p = a^p. 

Then, the equation of state p = Ke^'^~^^ ^p"', where s - entropy in the units of the universal gas 
constant 3fJ and squared sound velocity = yp/ p, gives the following differential relations: 

da^ , , 9 7 — 1 (is , , 
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1 dp Qeff -y — Ids 
pdr 0? 7 dr 



(IH]) shows that squared buoyancy frequency is simply 

7-1 
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9eff 



ds 

dr ' 



(9) 



So for stratified medium the eq. ^ is modified just by the additional term proportional to the 
entropy gradient. Finally, from the system ^ we obtain the equation for infinitesimal perturbations 
with respect to p: 



d 
dr 



m 



D-rP 



rp d I 
D dr 

(cu — mVtY j 



2m Id /Vtp 
uj — mVt \ dr \ D 



2VLp 



1 1 



.9eff, 



D Qeff J 

d (rp N"^ 



+ 



dr\D Qeff, 



^ = 0. 

p 



(10) 



where D = Dh + N'^, 



K 



{uj — mQ)"^ - the quantity signed as "D" for homoentropic flow 



in ZS2. The squared epicyclic frequency is the same as before: 

r dr ^ / 

Excluding the difference in notations, the main equation (|TOl) coincides with one used by Lovelace 
et al. (1999) for perturbations integrated in the z-direction. For A^^ = the eq. (|TOl) comes to the 
equation for the Eulerian perturbation of enthalpy in homoentropic fluid used in ZS2. In the limit 
7 ^ oo the equation ffTOj) determines the behavior of small perturbations in stratified incompressible 
flow. At the same time the entropy gradient is connected directly with the variable density: 



ds 
dr 



1 dp 
p dr'' 



and the squared buoyancy frequency is: 



9eff 



1 dp 
p dr^ 



:i2) 



At last, we write here the expressions for v^. and v^p, in terms of p/p: 

Vr = -p: (u; - mil) -r { - ] - { (^ - ™') " > 

D [ dr \pj \ r Qeff J \pj\ 

^ D[2n\dr\p) 9eff\p)) r\ lu - mil J \p J _ 

2.2 The boundary conditions 

We have to impose boundary condition for p/p. At the free boundary for perturbed as well as for 
unperturbed flow h = = p = 0. Here - is an enthalpy of perfect gas. This implies vanishing of 
the Lagrangian perturbation of enthalpy at the boundary of the unperturbed flow Ah\b = 0. Further, 
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for the fixed fluid particle dh = dp/p + TdS = dp/ p since the motion is isentropic. Consequently, 
A/i = /S.p/ p = at the free boundary. Using the relation between small Eulerian and Lagrangian 
quantities (look Tassoul, 1978) we have: 

^ + &i^l=0 (13) 

p p dr 

The eq. ( fT3ll yields the boundary condition for p/ p. Here it's necessary to use the relation between 
and Vr have been used to obtain ([6]) and expression for Vr which was written at the end of the last 
section: 




ri,r2 



(14) 



One can check that (HM is exactly the same as the boundary condition for the homoentropic flow 
used in ZS2. In particular, the eq. fll4p doesn't contain terms with variable entropy. 



3 Angular velocity and entropy profiles 

The rotation curve is determined by gravitation and the pressure gradient. In the previous research 
(ZSl and ZS2) we used the two profiles of angular velocity of the basic flow: 

HM^n. ((^)-%-^.„(2.i^))"^ (15, 

f!(r) = f!„ (-!-) ' (16) 

here ri and r2 - are the boundary points of the flow, ro - point where the matter rotates with the 
Keplerian frequency Qq- The coefficient K sets the deviation of Q{r) from the Keplerian profile. All 
frequencies will be in units of Qq and all the distances - will be given in units of tq below. To define 
the radial extent of the fiow we'll use the quantity w = {r2 — ri)/ro x 100%. 

The power law (fT6|) has been traditionally used in the research of axisymmetric flow instabilities. 
This profile is characterized by the maximum of Qejf at the edge of the fiow for q > 3/2. The rotation 
according to (1151) implies that Qeff obeys the sinusoidal law and vanishes at the boundaries. In ZS2 it 
was showed that Qeff — >■ in ri and r2 entails the absence of growing sound modes. For consistency 
we take here the same rotation profiles although the main attention will be paid to the power law. 

To solve the equation (ITOl) one should calculate first a^(r). Integrating the first expression in ([H]) 
we have: 

^2 ^ g{7-i)./7 / _ 1) e-(^-^)^/^ dr +C^ , (17) 



{l~l)geff e-(^-^)^/^ dr+C, 

where Qeff is defined by the expr. ([71), and Ci for specific w is determined from the condition that 
must vanish at ri and r2. Let's notice that for profile f|T5|) the boundary points are defined beforehand, 
since tq = (r2 + ri)/2. For (ITUI) ri and r2 must be calculated together with the constant Ci. 

For the certain parameters a^(r) was determined by the Chebyshev polynomials approximation of 
the tabular function obtained by the numerical integration of the equation fll7p . Note that we could 
avoid difficulties with calculation using an appropriate angular velocity profile. But in this case 
we would miss the possibility to compare the new results with the previous ones found in SZl and 
ZS2 in the limit of uniform entropy distribution. 
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Besides, we'd lose a chance to explore the effects induced by stratification purely, i. e. with the 
same distribution of effective gravity and vorticity in the basic flow. 

In the limit of 7 — oo one doesn't have to find explicit dependencies of the flow quantities (such 
as p or p) on r to solve ( |TOl) . However, in case of power law rotation profile we need to determine 
the boundary points, what can be done using the stationarity condition. Since ge/j = we have: 

p = J geffpdr + C2 (18) 

To integrate the equation (|T8i) . one should know p(r). For the incompressible stratified medium we 
get p(r) by integration of the equation (ITT]) with the certain entropy profile. 



Entropy profile 

The entropy distribution was defined as follows: 

s(r) = So - si (r/ro - l)^ (19) 

where as before, tq - is the point in the flow with Keplerian rotation, so that Qefj = 0, sq, si - are the 
constants. This form is convenient because when si 7^ the entropy gradient is directed in the same 
way with respect to ge// everywhere in (ri, r2), in the other words A^^ has the same sign going to 
the zero in the vicinity of tq. Moreover, the function (JT9l) is monotonous on the each side of tq, what 
allows to set entropy before calculation of the boundary points necessary for the power law rotation 
profile. Note that for the fixed w the result of calculations doesn't depend on sq. Indeed, adjusted 
for the dep. ffT^ squared sound speed ffTTj) is: 

a\r)=Fi{r) + L{so,Ci)F2{r) 

Now, if in the certain boundary point = 0, we find L{so, Ci). It means, that sq controls only Ci, 
and leaves to be unchanged. Then, since the equation ( fTOj) contains only the derivatives of s(r), 
we concluded that eigenvalues of the boundary problem are not affected by the value of Sq. So Si is 
the only parameter that characterizes the stratification of the flow. 



Opposite to non-axisymmetric an axisymmetric perturbations obey the Heiland's criterium (Tas- 
soul, 1978). Namely, for considered stationary flow the necessary and sufficient condition of it's 
stability to perturbations with axial symmetry is the inequality: 

K^ + N^ >0 (20) 

Thus the stratiflcation with Si < contributes to the stabilization of the flow (to axisymmetric 
perturbations), and vice versa, with si > - to the instability in comparison with the homoentropic 
case. For Q{r) oc r~'^ with q < 2 and Si > the flow with sufficiently small w will always satisfy to 
Heiland's criterium. However, the increase of w will cause the break of (1201) close to the boundaries, 
since the entropy gradient goes up while moving away from tq. 

One more criterium is valid for stratified fiows. It is the Richardson criterium. We've mentioned 
the authors who generalized it for baroclinic axisymmetric fiows. For barotropic configuration we 
consider here it yields that 

A^^ 1 

Ri = ^ > - (21) 

{rdn/drf 4 

everywhere in (ri, r2) is the sufficient condition of the stability (see also Glatzel, 1990). Ri - is the 
Richardson number. It's clear that (12 ip is never fulfilled here, since for any si i?i — >■ while r — tq. 
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4 Numerical calculation 



We'll solve the boundary problem consisting of the equation ffTOl) and boundary conditions f|T^ at the 
interval (ri, for compressible and incompressible fluid. Let's rewrite ffTUl) in the form convenient 
for integration: 




D' 


Ar2 




+ 


15 


9eff 


D' 




15 


9eff 



« \9effJ 

with the finite compressibility and 

p\ I i U' N"^ 



9eff 




+ 



(22) 





2m 



UJ 



m 



rriQ 




[UJ 







(23) 



\9effJ 9eff \r 

in the incompressible limit. In the equations (1221) and (1231) the derivative on r is marked by a stroke. 
Since we're interested in the growing modes the integration can be implemented along the real axis 
according to the Lin's rule (1945). Each of the equations (1^ and ([23D have been separated into 
the real and imaginary parts, so that the system of four first order equations consisting of the real 
quantities have been integrated. 

For Ui the equation (l23l) doesn't contain singular points at real axis so the boundary problem 
have been solved as in ZSl. The four specified linearly independent vectors have been used as initial 
conditions in ri to obtain the corresponding solutions in r2. The boundary conditions in ri and r2 
give then four algebraic equations and the eigenvalues should be obtained by setting the determinant 
of corresponding matrix equal to zero. The finite compressibility involves the sound speed to vanish 
at the boundaries. The main equation contains now singularities in ri and r2 (for power law rotation 
it's the first order poles and for ffT^ profile it's the second order poles) so the numerical algorithm 
have to be adjusted appropriately. 

As in ZS2 the solution should be expanded to the generalized series in the vicinity of ri and r2: 



p 
P 



ai{r - ri,2) 



(24) 



Substituting ( |2^ into (1221) and multiplying it by (r — ri 2)^ one gets the recurrent expressions for 



and the squared equation for /i (see ZS2). The last one defines two independent solutions of (l22l) 
close to the boundaries with some /xi and fi2- The regularity condition assigns what solution should 
be chosen. After that the numerical integration from left and right boundary gives two solutions in 
some point inside (ri, r2). The coupling condition can be rewritten as the determinant of the certain 
matrix equal to the zero as above and again gives an eigenvalues. At that fii and /i2 turn out to be 
independent on the entropy distribution in spite of the fact that the recurrent expressions for differ 
from that in ZS2. For the power law rotation /i = as before corresponds to the regular solution 
and the relation between oq (value of p/ p at the boundary) and oi (value of (p/ p)' at the boundary) 
as before coincides with the boundary condition ( |T4l) (As it's shown in ZS2 for the uniform entropy 
distribution). 
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For rotation with quasi-sine deviation from the Keplerian law we choose again the solution with 
Re{fi) > 1, what simultaneously implies that it satisfies to the boundary condition when 5'e//|ri,2 = 0. 

It's appropriate to mention here one popular way when the stratification of basic flow is set by 
the relation p oc p^, where F 7^ 7 (see, for instance, the paper by Kojima, 1989). This approach 
can't be implemented here since at the boundaries r — > ri^2, P 0, and entropy and its derivatives 
are not flnite in ri^2- So the coefficients of f l22l) have the singularities also because of the terms that 
contain entropy gradients. Then it's easy to demonstrate that, at least for the power law rotation, 
the relation between and ai with /i = won't be equivalent for (fT4|) . so that the boundary problem 
doesn't have the solution. The boundary condition is equivalent to the regularity condition only if 
s(r) with it's derivatives is finite at the boundary. 

5 The results 

5.1 Instability of flow with the power law rotation profile. 

There are growing sound and surface gravity modes in axisymmetric fiow with free boundaries (Pa- 
paloisou & Pringle (1984, 1985, 1987), Blaes & Glatzel (1986), Goldreigh et al. (1986), Glatzel 
(1987a,b), Sekyia & Miyama (1988), Kojima (1989) and others). An entropy distribution is sup- 
posed to modify mentioned branches of instability and possibly the growing internal gravity modes 
should appear (Glatzel, 1990, 1991). The results of calculations will be presented here mostly at the 
plots, where the curves of increments will be drawn. We also won't limit ourselves with the region 
stable to axisymmetric perturbations according to the Heiland's criterium. In the curves of uJi{si) 
the transition from one region to another will be marked by a stroke. 

Calculations with finite compressibility In fig. 1 uJi{si) and UJr{si) are displayed for the growing 
surface gravity mode with m = 10, which form so called principal or incompressible branch of 
instability, since uoi > exist also in the limit 7 — > 00. In the absence of stratification for q < 2.0 
and 7 < cxo an increment of principal branch disappears at some maximum value of w. As we can 
see for small values of w an increment monotonously increases while si gets higher, i.e. while the 
configuration becomes less stable according to the Heiland's criterium fl20|) . However, for greater w 
minimum appears in the increment curve and then the growing incompressible mode disappears at 
certain si > 0. For even greater w surface gravity mode stops growing in the fiow with s = const, but 
at the same time increment emerges both for the positive (si < 0), and the negative (si > 0) entropy 
gradient. The last result is consistent with the calculations of Kojima et al. (1989), who made the 
conclusion that entropy growth in the opposite direction to effective gravity (here it's means si < 0) 
contributes to widening of instability range in the values of w. The similar effect gives the calculation 
of uji{si) for m = 1 and different values of q (look at fig. 2). One can see that while the vorticity 
gradient is increasing (i.e. q is decreasing) stabilization sets first in the homoentropic case. Beside 
that, in fig. 2 dashed lines denote the additional branches of instability. As it'll be clear below, these 
are growing internal gravity modes caused by stratification. However, these modes were found in the 
present work only in the region unstable with respect to axisymmetric perturbations. 

In a more intricated manner stratification affects the growing sound modes. All the branches of 
sound modes increments for specific w are depicted in fig. 3. We've used here the results of ZS2, 
where the detailed calculation of sonic instability was implemented, exactly, the calculation oiLJi{w), 
with m = 10 and q = 1.58. We've used the value of w corresponding to the modes coupling for one 
of the increments. In the homoentropic fiow the modes coupling entails the strong increase of Ui in a 
narrow range of w. In fig. 3 the modes coupling with si = occurs for branch (5). To the right side 
from the vertical dashed line the fiow is unstable according to the Heiland's criterium. The increment 
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Figure 1: Increment (pic. (a)) and pattern speed (pic. (b)) of non-axisymmetric modes of perturbations as a function 
of entropy gradient for the power law rotation profile. The numbers are the values of w in percents. The solid lines 
denote the growing surface gravity modes, the dashed lines - the sound modes. The strokes in the curves mark the 
edge of stability to axisymmetric perturbations. For curves corresponding to w = 9%, 10%, 15%, with si > the flow 
is unstable to axisymmetric perturbations. 



10 



(a) 



CO ; 



0.5 



0.4 



0.3 




(b) 



0.95 



CD, 




0.85 



Figure 2: The same as in fig. 1, but for fixed w and the dashed lines denote now the growing internal gravity modes. 
The numbers arc the values of q. For curves with q = 1.6, 1.7 the flow in unstable to axisymmetric perturbations. The 
curves of pattern speed of internal modes are not depicted because > 1- 
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branches fall into three categories: (1),(8),(2),(7),(3) - the perturbations growth disappears both at 
some si < and at some si > 0; (4), (6) - increment gradually increases with si, (5) - corresponds 
to the modes coupling and uJi{si) behaves in the most tangled way. Exactly, at some value of 
Si = sf > increment (5) vanishes, but comes back at once and grows then with si. The same 
features of modes coupling arise in dependence uji{w) for homoentropic flow presented in ZS2. To 
check if the presence of s = sl^, where cUj — is the common peculiarity for modes coupling, we've 
calculated few a;j(si) for different w. In can be seen that all determined uji{si), have a break point 
in increment curve (fig. 4). 

At last, branches of sonic instability for m = 1 one can see in fig. 5. The calculations here were 
mostly carried out in area unstable to axisymmetric perturbations. It was revealed that both at 
some Si < and at some Si > increment vanishes. In ZS2 we discussed that the growing sound 
modes with m = 1 arise as a result of two mechanisms: resonant amplification by the basic flow in 
critical layer (mechanism Landau) and coupling with decaying surface gravity mode, what sets it 
apart from the above instability being the consequence of coupling of two sound modes. Clearly, the 
increment behavior differs from what we've seen in fig. 4. The main is that in this case for none of 
the branches increment arises again after vanishing at some Si > s^^. 

While calculating the presented curves we tried to find the growing internal gravity modes in area 
stable to axisymmetric perturbations. However we failed to do that and the only one instability that 
arises due to stratification is denoted by dashed lines in fig. 2. We've made an attempt then to find 
growing internal gravity modes in the incompressible limit, since this branch of instability must be 
insensitive to the compressibility. 

Calculations in incompressible limit An incompressible limit was used here first of all as an additional 
way to check the results presented so far, since the numerical algorithm differs considerably in that 
case (look above). In fig.6 we put increment curves with dependence on 7 (pics a & b) and with 
dependence on si ( pics c & d) calculated in approximation concerned here. For 7 = 5/3 the 
corresponding values of Ui and 00^ can be found in fig. 2. The decrease of sound speed suppresses an 
internal gravity modes growth (dashed curves) and its increment vanishes at some 7 > 1. On the 
contrary, the surface gravity modes are unstable right up to 7 — > 1. Then, as it must be, while Si ^ 
uji of surface mode diminishes up to it's value in homogenerous incompressible fluid (Jaroszynski, 
1988, SZl). Concerning internal mode it's growth disappears at some si 7^ 0. 

The parameters in fig.6 correspond to unstable axisymmetric perturbations. We've already noticed 
that we missed to find growing internal gravity modes in the flow with free boundaries and 7 < 
00 that is stable according to Heiland's criterium. The numerical analysis confirmed this also in 
incompressible limit. However, it turned out that if someone impose rigid boundaries the increments 
of modes that grow with free boundaries increase significantly exceeding Qq. 

Moreover, rigid boundaries involve the instability of configuration stable according to Heiland's 
criterium fl^Uj) . Note that Glatzel (1990) also have found growing internal modes in a problem 
with rigid boundaries. At a glance, such influence of boundaries on the unstable internal modes is 
embarrassing and the stratification must play the main role. However one should remember that 
we consider the amplification of internal oscillations by the shear flow. The neutral internal modes 
certainly exist independently on the kind of the boundaries. But the ability to grow (to decay) for 
example owing to interaction with the basic flow is controlled not only by a vorticity gradient at 
the critical layer (as it's for homoentropic flow) but additionally by the shape of perturbations field 
itself as it was shown, in particular, by Troitskaya & Fabrikant (1989) (see also the monograph by 
Stepanyants & Fabrikant, 1996). The latter is confirmed also by Glatzel's paper (1990), where he 
found that the internal modes grow (decay) outside the regions of coupling what can happen only 
if perturbations resonantly interact with the basic flow. But he used the simplified angular velocity 
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Figure 3: The same as in the fig. 1, but for fixed w and the increment and pattern speed of sound modes only are 
presented. The numerals are the numbers of instability branches. The vertical line marks the limit of stability to 
axisymmetric perturbations. 
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(a) 





Figure 4: The same as in the fig. 3 , but for different w, signed each curve. The branches presented are the 

result of sound modes coupling (for si = 0). The strokes in the curves denote the edge of stability to axisymmetric 
perturbations. For w = 33.15%, 27.7% the flow is stable to axisymmetric perturbations in the whole range of si. 
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Figure 5: The same as in the fig. 3 but for m = 1. 
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Figure 6: In pics, (a) and (b) - the dependences cui{^) and coril), in pics. (c) and (d) - the dependences LOi{si) and 
CiJr{si). The solid Hnes denote the growing surface gravity modes, the dashed hnes denote the growing internal gravity 
modes. 
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Figure 7: Increment of the surface gravity mode as a function of entropy gradient for the Keplerian rotation with 
quasi-sine deviation. The angular velocity profile corresponds to the profile of ^l{r) in the case of power law rotation 
with the same w (look comments in text). In the marked range of si the flow is stable to axisymmetric perturbations. 

profile corresponding to the constant vorticity in the flow, what excludes the interaction of small 
perturbations with the flow. 

The perturbations proflle which is responsible for its growth along with the vorticity gradient in 
stratifled fluid is the solution of the boundary problem, so that the ampliflcation of internal modes 
is determined also by the boundary conditions. 

5.2 Instability of rotation with quasi-sine deviation from the Keplerian profile 

The instability of flow with such proflle was investigated in the previous papers ZSl and ZS2. Par- 
ticularly, it was revealed that the stabilization always occurs at some least i^' > (look the formula 
( IT5l) ), and for j < oo there's no sonic instability. The last implies that in non-stratifled flow only 
surface gravity modes can grow. Accordingly, in the flg.7 we can see how this sort of perturbations is 
modifled by the stratiflcation. The angular velocity proflle is set in order to make equal the maximum 
enthalpy in the flow with si = to maximum enthalpy in the flow with the same radial size w but 
rotating with a power law proflle with q = 1.58 (see ZS2). Note that calculations were carried out for 
7 = 2.5, since the flow turns out to be stable for ordinary 7 = 5/3. Several curves uji{si) correspond 
to the different values of w. The dependence on si is trivial: uJi increases while the flow becomes 
less stable according to the Heiland's criterium ( |20|) . Similar to the power law rotation the negative 
entropy gradient (si > 0) widens the range of w when the growing surface gravity modes exist. 

6 Conclusions 

Well, we investigated in the two-dimensional approximation the instability of the laminar axisym- 
metric stratifled flow with free boundaries to inflnitesimal non-axisymmetric perturbations. In the 
case of the power law rotation proflle the essential influence of entropy gradient on the growing sound 
and surface gravity modes was revealed. With both negative and positive sign of entropy gradient 
relative to ge// direction the range of w and q where growing surface gravity modes exist becomes 




17 



wider. Specifically, its increment vanishes at greater w and smaller q. For sound modes the incre- 
ments fall into three categories. In the first one Ui gradually increases while the flow becomes less 
stable according to the Heiland's critcrium (i.c while Si gets higher). In the second case increment 
vanishes both at some positive (si < 0) and negative (si > 0) entropy gradient. At last, the third 
category contains increments being the result of modes coupling. Corresponding u!i{si) has a break 
in s1^, where oOi — > 0. Similar peculiarities in increment behavior has been revealed also in ZS2 for 
u}i{w) which also had the points cjj 0. 

Besides the sound and surface gravity modes we tried to find growing internal gravity modes. The 
latter was done as well in incompressible limit. It turned out that mentioned sort of perturbations 
exists only in the flow that is unstable according to the Heiland's criterium, i.e. unstable to axisym- 
metric perturbations. However, with rigid boundaries the growing internal modes where found also 
in the stable to axisymmetric perturbations region (see also Glatzel, 1990). We suppose such crucial 
role of the boundaries to be the consequence of specific character of amphfication in a stratified 
medium. In particular, the transition of energy to perturbations from the basic flow is determined 
now not only by a vorticity profile in the critical layer but also by the profile of perturbations, which 
in turn depends on the boundaries (Troitskaya & Fabrikant, 1989). 

In the case of the Keplerian rotation with quasi-sine deviation we found that the increment of the 
surface gravity mode increases gradually while si gets higher. 

We also discussed in the paper that it's incorrect to set the stratification in a problem with free 
boundaries by a polytropic law with the index different from the adiabatic value, since in this case the 
perturbations don't satisfy the boundary conditions. The last is the consequence of infinite entropy 
and it's derivatives in the boundary points. It reveals that the condition of perturbations regularity 
at the boundaries is equivalent to the boundary condition (i.e. Ap — 0) only if entropy and it's 
derivatives have finite values at the boundaries. 

This paper was supported by grant RFFI-NNIO 06-02-16025. 
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